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Optimal Control of Multiphoton Ionization Processes in Aligned I 2 Molecules 
with Time-Dependent Polarization Pulses 
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Multiphoton ionization processes in aligned I 2 molecules are actively controlled by the homemade 
pulse shaping system, with which a time-dependent polarization pulse can be generated and controlled 
We find a correlation between a femtosecond time-dependent polarization pulse and the production 
efficiency of evenly or oddly charged molecular ions. We achieve much better controllability of the 
correlation with a time-dependent polarization pulse than with a pulse having a fixed ellipticity. The 
results suggest the existence of an unknown tunnel ionization mechanism which is characteristic of a 
time-dependent polarization pulse. Our experiments point to new directions in optimal control studies 
with molecular systems, as discussed in the text 
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Based on the combination of the femtosecond pulse 
shaping technique [1] associated with the state-of-the- 
art ultrafast laser technologies and optimization algo- 
rithm, it has become possible to optimize quantum 
processes in atoms and molecules and achieve the desired 
results with a learning-loop pulse shaping system [2], 
This approach is versatile and applicable even to the 
laser-molecule system for which the accurate Hamil- 
tonian is not known. It corresponds to the experimental 
solution of a time-dependent Schrodinger equation by the 
use of measured values and the subsequent feedback with 
an optimization algorithm. The meanings of such a class 
of study are summarized as follows: (1) to elucidate the 
interaction between an ultrafast laser field and a mole- 
cule by analyzing the optimized pulse, (2) to verify 
hypotheses, (3) to search new phenomena by optimizing 
a specific experimental parameter, and (4) to optimize 
experimental conditions. These meanings are comple- 
mentary with one another. 

Since its first proposal of feedback quantum control [2], 
many demonstrative experiments have been performed 
for molecular systems (Refs. [3-5], and references 
therein). However, in experiments performed in gas- 
phase molecules, they are randomly oriented and give 
only the averaged results concerning the molecular ori- 
entation. This limits the outcome of optimal control ex- 
periments and prevents us from studying the correlation 
between the laser polarization and the molecular axis. In 
fact, there are many important polarization dependent 
processes such as parallel (perpendicular) transitions 
and enhanced ionization [6-8]. Here we demonstrate 
optimal control experiment where a sample of aligned 
molecules [9,10] is employed for the first time to our 
knowledge. Another important feature of the present 
experiment is that a time-dependent polarization pulse 
[11,12] is used as a new control parameter. Better con- 
trollability obtained with a time-dependent polarization 
pulse than with a pulse having a fixed ellipticity suggests 



the existence of an unknown underlying physics and can 
trigger theoretical investigations to elucidate it Actually, 
this is the case of the present study. A time-dependent 
polarization pulse should be a powerful tool especially for 
a sample of aligned molecules. A special emphasis is 
placed on the fact that the present experiment is also the 
first demonstration where both external (rotation) and 
internal (ionization and excitation dynamics) degrees of 
freedom in molecules are simultaneously controlled. Our 
experiment should point to new directions in feedback 
quantum control and molecular physics experiments as 
discussed below. 

Here we pay close attention to the multiphoton ion- 
ization processes in I 2 molecules. They have been 
theoretically and experimentally studied over decades. 
The molecule can ionize sequentially [13] or nonsequen- 
tially [14-16]. The nonsequential double ionization is 
caused by recollision: (1) First an electron tunnels 
through the potential modified by the strong linearly 
polarized laser field, (2) then the electron is driven by 
the field and returns to the parent ion after the field 
reverses its direction, and, finally, (3) in the electron- 
ion recollision that might follow, one of the remaining 
electrons is collisionally ionized if the recollision elec- 
tron has a sufficient kinetic energy. If the polarization 
deviates very much from linear, when the electron re- 
turns to the ion, it is offset in the lateral direction and it 
can miss. As a corollary of nonsequential double ioniza- 
tion, we can expect that it contributes more to the pro- 
duction of evenly charged molecular ions such as I 2 + than 
to that of oddly charged molecular ions such as I 2 and 
I^ + . It means that we can expect a correlation between the 
polarization of the laser field and the production effi- 
ciency of evenly or oddly charged molecular ions. 

We experimentally investigate the correlation. Here an 
outline of the experimental setup (Fig. 1) is described and 
a fuller description can be found in our recent papers 
[17,18]. A pulsed supersonic beam of I 2 molecules is 
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introduced into a time-of-flight (TOF) mass spectrometer 
by expanding I 2 molecules (~ 1 Torr) buffered with argon 
gas (760 Torr) through a 0.25-mm-diameter nozzle. The 
molecular beam is parallel to the TOF axis and crossed at 
90° by the focused laser beams. In order to align mole- 
cules, we use pulses from an injection-seeded Nd:YAG 
(Nd-doped yttrium aluminum garnet) laser (Quanta-Ray, 
GCR-130, wavelength A = 1064 nm, pulse width r~ 
12 ns, and the peak intensity / ~ 2 X 10 12 W/cm 2 ). The 
polarization direction is set to be parallel to the TOF axis 
and thereby the I 2 molecules are aligned along the TOF 
axis. Under the present conditions, the ensemble average 
of the squared alignment cosine is estimated to be —0.7 
by the two-dimensional ion imaging technique [9,10]. 
Intense femtosecond Tiisapphire laser pulses (Spectra- 
Physics, Super Spitfire, central wavelength A ~ 800 nm, 
r ~ 45 fs, and / — 2 X 10 14 W/cm 2 ) are used to ionize 
the I 2 molecules. We carefully adjust the focal spot size of 
the femtosecond pulse to be smaller than that of the YAG 
pulse to ensure that we probe only aligned molecules that 
have been exposed to the YAG pulse. The fragment ions 
produced from photodissociation and Coulomb explosion 
are accelerated by a static electric field toward a micro- 
channel plate (MCP, Galileo, AP TOF- 18) detector posi- 
tioned on axis with the TOF axis. 

We adopt a fitness function F defined by F = /(l, 1)/ 
[1(1 0) 4- I(i t 2) + 1(2, 1)], where /(m, n) stands for the 
integrated I m+ signals of the fragmentation channel 
jm+ + jn+ produced from l^ +n)+ . In general, higher laser 
intensities are needed to produce more highly charged 
molecular ions. In order to reduce such an intensity- 
dependent effect on the fitness F, we compare ion signals 
from both I 2 and I 2 + (oddly charged molecular ions) with 
those from I 2 + (evenly charged molecular ions). With a 
sample of aligned molecules, two fragmentation channels 
(1,1) and (2,0) produced from I^ + should be treated 
separately as discussed below. 

The peak intensity of 2 X 10 14 W/cm 2 (linear polar- 
ization) gives the ponderomotive energy of U p = 12 eV. 
Depending on the phase of birth of the electron, it can 
recollide with its parent ion with the maximum kinetic 
energy of 3AlU p = 38 eV [14,15], which is enough to 
knock the remaining electron from I2 (ionization poten- 
tial —20 eV) to produce I 2+ by nonsequential double 
ionization caused by recollision. Using such femtosecond 
pulses, we measure the ellipticity dependence of the fit- 
ness F with the combination of a quarter-wave plate and a 
half-wave plate. The major axis of the elliptical polariza- 
tion is kept parallel to the TOF axis, i.e., the molecular 
axis of aligned I 2 molecules. (See Fig. 1 and the orienta- 
tion angle 6 = 0 in the present case. The orientation angle 
6 is defined by the angle between the major axis and the 
TOF axis.) This is the basis of including (1, 1) and ex- 
cluding (2, 0) in the fitness F as the fragmentation channel 
from I^ + . We are interested in the contribution of the 
nonsequential double ionization to the production of 
l2 + , and the (1, 1) channel is expected to be dominant 
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FIG. 1 (color online). An experimental setup. A pulse shaper 
is used to generate tailored femtosecond pulses for learning- 
loop optimal quantum control. Experimental feedback signals 
are collected from time-of-flight (TOF) mass spectra that 
result from the interaction of the tailored femtosecond pulses 
with aligned molecules. These signals are processed in a com- 
puter to iteratively improve the applied pulse shape based on 
genetic algorithm until an optimal result is obtained Ellipticity 
e and the orientation angle 0 are also illustrated 

when the laser polarization is parallel to the molecular 
axis. This might be explained by the following simple 
consideration. An electron initially localized in one of the 
I atoms tunnels through the lower and thinner potential in 
between the two I atoms. The electron is driven by the 
laser field and returns to and recollides with the parent 
ion after the laser field changes its direction. Then the 
electron localized in another I atom should be more likely 
to be ionized because the slope of the molecular potential 
combined with the laser field is also reversed In this case, 
the resultant fragmentation channel should be (1, 1). On 
the other hand, when the laser polarization is perpendicu- 
lar to the molecular axis, the ionized electron should be 
more likely to recollide with the I atom from which the 
electron was ionized In this case the resultant fragmen- 
tation channel should be (2, 0) [19]. 

The ellipticity e is defined by the ratio of the minor 
axis to the major axis of an elliptical polarization and 
varies from e = 0 (linear polarization) to s = 1 (circular 
polarization). We accumulate ion signals with 2000 laser 
shots for each ellipticity and the results are presented in 
Fig. 2. The fitness F is 0.215 for linear polarization and 
gradually decreases down to 0.195 for circular polariza- 
tion. It means that the relative production efficiency of 
evenly charged molecular ions (I^) is high for linear 
polarization and that of oddly charged molecular ions (I2 
and I^ + ) is high for circular polarization, which is quali- 
tatively consistent with our hypothesis. The contrast de- 
fined by (the maximum fitness)/(the minimum fitness) is 
—1.1. When we employ the (2, 0) channel instead of (1, 1) 
in the fitness F it presents opposite ellipticity dependence 
and increases as the ellipticity e is increased. This obser- 
vation can be explained by considering that the time 
period during which the electric field vector of the pulse 
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FIG. 2. Ellipticity dependence of the fitness F defined 

in the text The contrast defined by (the maximum F)/ 
(the minimum F) is —1.1. 



is perpendicular to the molecular axis contributes to the 
production of the (2, 0) channel. These observations al- 
ready demonstrate the efficacy of a sample of aligned 
molecules. When we include both the (1, 1) channel and 
the (2, 0) in the fitness F, we cannot observe any clear 
ellipticity dependence. This may be due to the mutual 
compensation of two different situations discussed in the 
previous paragraph, which can be validated by all these 
observations. 

Then we optimize the fitness F with the homemade 
learning-loop optimal control system where a time- 
dependent polarization pulse can be generated and 
controlled [11,12]. We use a typical 4/ configuration 
pulse shaper but remove polarizers to generate a time- 
dependent polarization pulse. We employ the genetic al- 
gorithm (GA) [20] as an optimization algorithm and 
search for an optimal time-dependent polarization 
pulse. We use the ellipticity for each frequency component 
as the only parameter in GA, and the phase for each 
frequency component is kept constant As an initial pulse 
used in GA, we use a "randomly polarized pulse" for 
which the ellipticity for each frequency component is 
randomly defined With a sample of aligned molecules, 
fragmentation signals in the TOF spectrum become 
sharper compared to those with randomly oriented mole- 
cules. Sharper and isolated signals are very advantageous 
to be used in the fitness F for feedback learning loop. 

We maximize and minimize the fitness F and its 
evolutions are shown in Fig. 3 as a function of generation 
of GA. The maximum (minimum) fitness obtained is 0.34 
(0.14) and the contrast is —2.5, which is much larger than 
that obtained with a fixed ellipticity (1.1). Our results 
show a much better ability to control the charged state 
of molecular ions with time-dependent polarization 
pulses than with elliptically polarized pulses having a 
fixed ellipticity. We characterize the optimal pulses with 
the technique known as POLLIWOG (polarization 
labeled by interference versus wavelength of only a glint) 
[21]. The results are shown in Fig. 4. The intensity profiles 
of the optimized pulses are confirmed to be almost iden- 
tical with those used in Fig. 2. The polarization is pre- 
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FIG. 3. The fitness F defined in the text is maximized and 
minimized as a function of generation of the genetic algorithm 
The contrast defined by (the maximum F)/(the minimum F) is 
—2.5, which is much larger than that obtained with a fixed 
ellipticity as shown in Fig. 2. 

scribed by both the ellipticity e and by the orientation 
angle 0. When the fitness F is maximized, i.e., the pro- 
duction efficiency of l\ + is optimized with that of oddly 
charged molecular ions suppressed, the pulse becomes 
linear polarization parallel to the TOF axis (0 = 0) at 
around the peak of the pulse. Such a polarization corre- 
sponds to the one which we consider to be a favorable 
situation for the production of the (1, 1) channel due to 
recollision. The peak intensity of 1.3 X 10 14 W/cm 2 
gives the maximum recollision energy of 3AlU p = 
24 eV, which is enough to knock the remaining electron 
from I2 to produce I 2 ^ by nonsequential double ioniza- 
tion caused by recollision. On the other hand, when the 
fitness F is minimized, i.e., the production efficiency of 
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FIG. 4. The time-dependent polarizations and the intensity 
profiles of the optimized pulses characterized by the homemade 
POLLIWOG system In the upper panels, both ellipticity e 
(solid curve, left ordinate) and the orientation angle 0 (dashed 
curve, right ordinate) are presented. Ellipticity e and the 
orientation angle $ are illustrated in the inset The sign of e 
is related to the helicity of the elliptical polarization. If positive 
values of e define left elliptical polarization, negative values of 
e define right elliptical polarization and vice versa The lower 
panels show the instantaneous intensity profiles expressed as 
/W = /o/W> where / 0 and f(t) are the peak intensity and the 
pulse time profile, respectively. The left and right panels 
correspond to the optimized pulses which give the maximum 
and the minimum F, respectively. 
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oddly charged molecular ions is optimized with that of 
l| + suppressed, the ellipticity stays small (e ~ 0.1) from 
-60 to +30 fs, i.e., over the main part of the pulse. Here 
we note that the orientation angle changes more rapidly 
than in the maximized case and it goes through zero just 
as the case with the fitness F maximized Qualitative 
features of the optimized polarizations are consistent 
with our expectations and the results shown in Fig. 2. In 
fact, the polarization becomes linear (elliptical) at around 
the peak (over the main part) of the pulse optimized to 
maximize (minimize) the fitness F. However, the ellip- 
ticity of ~~0. 1 over the main part of the pulse optimized to 
minimize the fitness F is rather small compared to the 
results shown in Fig. 2. The results suggest the impor- 
tance of the time evolution of the potential energy surface 
associated with the increase of the internuclear distance 
as well as the electrons' dynamics discussed above. 

Apparently, further theoretical studies are required to 
understand the underlying physics in the optimal control 
of the charged state of molecular ions. A tentative sce- 
nario is as follows: We note that the optimized pulses have 
polarization components perpendicular to the aligned 
molecular axis before (and after) the peaks in both cases 
of maximization and minimization of the fitness F. The 
polarization components perpendicular to the aligned 
molecular axis are preferable to a perpendicular transi- 
tion such as 4 3 Il lH [22]. When the molecule is 
excited on a dissociative potential such as A state and the 
internuclear distance increases to the critical region of 
enhanced ionization, the molecule is in a favorable situ- 
ation both for sequential and for nonsequential double 
ionization. In such a situation, there might exist a new 
ionization mechanism by which the ionization is en- 
hanced even with a relatively small ellipticity. The gen- 
eralization of enhanced ionization with a time-dependent 
polarization pulse, where electron localization, time evo- 
lution of internuclear distance, destructive interference 
between different excitation paths and/or molecular orbi- 
tals, and Landau-Zener transition should play important 
roles, is now in progress in our group and will be pre- 
sented elsewhere in the near future [23]. 

In conclusion, we find a correlation between a femto- 
second time-dependent polarization pulse and the pro- 
duction efficiency of evenly or oddly charged molecular 
ions using a sample of aligned molecules. Much better 
controllability of the correlation is achieved with a time- 
dependent polarization pulse than with a pulse having a 
fixed ellipticity. The novelty of our study is multifold and 
points to the following new directions in optimal control 
studies with molecular systems: (1) A sample of aligned 
molecules is used in the optimal control experiment for 
the first time. The use of aligned molecules enables us to 
find a correlation between the polarization of the laser 
field and the production efficiency of evenly or oddly 
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charged molecular ions, which should have been washed 
out with randomly oriented molecules. (2) In order to 
optimize the quantum processes in molecules, the time- 
dependent polarization pulse is applied to the learning- 
loop optimal control system for the first time. It does give 
a better controllability of the charged state of molecular 
ions than a pulse with a fixed ellipticity. The importance 
of the time-dependent polarization pulse is also pointed 
out to selectively prepare one of the enantiomers from the 
50:50 racemates [24]. (3) Using a sample of aligned 
molecules and a time-dependent polarization pulse, both 
external and internal degrees of freedom in molecules are 
simultaneously controlled 
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